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A brief review is given of some of the recent developments in the theoretical analyses 
of supersymmetric dark matter. These include the effects of uncertainties in the 
wimp velocity and wimp density and of the effects of uncertainties in the quark 
densities of the proton. Also analyzed are the effects of non-universalities in the 
gaugino sector and their effects on determining the nature of cold dark matter, i.e., 
if the neutralino is bino like, higgsino like, or wino like. The maximum and the 
minimum elastic neutralino proton cross sections are discussed and a comparison 
of the direct and the indirect detection arising from the capture and annihilation 
of neutralinos in the core of the earth and the sun is given. Some of the other 
recent developments are summarized. 



1 Introduction 

Because of the.^m:ent significant activity in dark matter searches on the ex- 
perimental sidcl'aQ there is renewed interest in the theoretical analyses of dark 
matter which are significantly more refined than in the previous years. Among 
the refinements is the inclusions of the effects of uncertainties in the input 
parameters in the theoretical predictions of event rates and of the neutralino 
proton cross sections as well inclusion of the effects of non-universalities, CP 
violating effects and the effects of coannihilation. The content of the paper is 
as follows: In Sec. 2 we discuss the effects of uncertainties in the analyses of 
dark matter. These include the effects of uncertainties in the wimp velocity 
and in the wimp relic density, and the effects of uncertainties in the quark 
densities in the proton in the analyses of dark matter. In Sec. 3 we give a 
discussion of the maximum and the minimum elastic neutralino-proton cross- 
sections. In Sec. 4 we discuss the effects of non-universalities and specifically 
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Figure 1: An exhibition of the effect of the variation of the wimp velocity and wimp relic 
density on the event rates in a Ge detector as a function of the fine tuning parameter <I>. 
Taken from Ref.^ 

the non-universalities in the gaugino sector on dark matter analyses. A brief 
discussion of the effects of /x on the composition of the neutralino and its role 
in determining the nature of cold dark matter, i.e., if it is dominantly a bino, a 
wino or a higgsino is given in Sec. 5. A comparison of the direct and the indirect 
detection of dark matter is given in Sec. 6. In Sec. 7 we give a discussion of the 
annual modulation effect in the direct detection of dark matter. In Sec. 8 we 
give a brief discussion of the effects of CP phases on dark matter. Conclusions 
are given in Sec. 9. 

2 Uncertainties in Theoretical Analyses of Dark Matter 

The direct detection of dark matter has been investigated by many authors ( 
see RefQ for some of the recent works). Several types of uncertainties enter in 
these analyses such as the effects of variations in the local wimp density, the 
effects of variations in the wimp velocity range, and the effects of uncertain- 
ties in the quark densities in the nucleons. Other effects not considered here 
are the uncertainties in the nuclear form factors, and effects of halo modelfl 
on the event rates in dark matter detection, p. We begin wiltL^a discussion of 
the analyses of uncertainties in wimp densitjQ and velocitju'ErEi. The current 
range of local wimp density lies in the rangctJ (0.2 — 0.7)GeVcm^^ . Defining 
^ = p^a / po one can parameterize the local wimp density in term of which 
for Po = 0.3GeVcm~^ gives 0.7 < ^ < 2.3. For the wimp velocity one typically 
assumes a Maxwellian velocilv distribution for the wimps. Estimates of the 
rms wimp velocity range givJEI v — 270 ± 50 km/s. The analysis is capied out 
within the framework of the minimal supergravity model (SUGRAp3. The 
soft SUSY breaking sector in the minimal version mSUGRA of this theory 
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is given by mg, r7ix/2i ^Oi and tan/? where, mg is the universal scalar mass, 
rfiii2 is the universal gaugino mass, is the universal trilinear coupling, and 
tan/3 =< H2 > / < Hi > where H2 gives mass to the up quarks and Hi 
gives mass to the down quarks and the leptons. The Higgs mixing parameter 
/I is determined by the constraint of radiative breaking of the electro-weak 
symmetry. It is also useful to define a fine tuning parameter $ so thatJ 
$ = (/i^/Mf + l/4)^/'^. The fine tuning parameter defines how heavy the 
SUSY spectrum gets. The effects of the uncertainly in the event rate as a 
function of the fine tuning parameter due to variations in the wimp relic den- 
sity and in the wimp velocity are shown in Fig.Jj. One finds that the effects 
of this type of uncertainty can lead to a variation in the rates by a factor of 2-3. 

Next we discuss the uncertainties in neutralipOiOmton cross-section arising 
from errors in the quark densities in the protoit3'ta'll3. The basic interaction 
governing the x~P scattering with CP conservation is the four Fermi interac- 
tion given by £e// = Xlt^laXQl^'i^PL + BPaJq + Cxxrriqqq + DxibX'>T^qQlbQ- 
For heavy target materials the neutralino-nucleus scattering is dominated by 
the scalar interactions which is controlled by the scalar x ~ P cross-section 
where 
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Here /^^ is the reduced mass and ff (i=u,d,s quarks) are the (u,d,s) quark 
densities defined by rripff =< p\mqiqiqi\p > (i=u,d,s). In the above C is 
the strength of the scalar interaction and consists of s channel contributions 
from the higgs h° , ° exchange and t channel contributions from the sfermion 
exchange so that C ~ C^o + Cho +C^. The uncertainly in the theoretical pre- 
dictions of ay-p{scalar) is dominated by the uncertainty in the quark densities 
ff. To study the uncertainties in ff it is best first to solve the quark densities 
analytically in terms of some judiciously chosen parameters. One find£j'EZI 

/r„.) = ^^(i±o^, = (2) 

iu,a; ^ _|_ 

where the parameters ^, x and CTjrJV are defined by 

< p\uu — dd\p > < p\uu + dd — 2ss\p > 



< p\uu + dd\p > ' < p\uu + dd\p > 

< p\2^^{mu + md){uu + dd\p >— (Tt^n (3) 
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Figure 2: The effect on ctx-p of the variations in the quark densities around the central 
value given in Sec. 2 



Similarly one can analytically solve for the quark densities in the neutron and 
the analytic relations provide an interesting connection between the quark 
densities in the proton and ia. the neutron. One finds that independent of the 



h&mfUl - fufd- One can use the 
to determine the ratio f/x. Om 



details of any input oi: 
baryon mass splittings 

Using the determination of x from lattice gauge analyse 
0.132 ± 0.035. Additional uncertainties can 
Here results from chiral perturbation theorjl 



alysis on the 
i/x = 0.196. 
one finds ^ = 
ifrom the quark mass ratios, 
give ^ = 0.553 ± 0.043, 



= 18.9 ± 0.8. Using the inputs above one finds 



0.021 ± 0.004, fl 



0.029 ± 0.006, /f = 0.21 ± 0.12 and = 0.016 ± 0.003, = 0.037 ± 0.007, 
/J* = 0.21 ±0.12. 



3 Maximum and Minimum Neutralino-proton Cross Sections 

We discuss now the numerical results of the ncutralino proton cross-sections 
with the quark densities as discussed in Sec. 2. In Fig. 2 cr^-p is plotted ex- 
hibiting the effects of the variations of the quark density as a function of m~^. 
The above analysis shows that the x~ P cross section cannot be computed to 
an accuracy of better than a factor of 3-5 with the current level of uncertainties 
in the input data. Fig. 2 also gives a plot of the maximum and the minimum 
of function of the neutralino mass. In the analysis of Fig. 2 we have 

allowed mg and r7ii/2 to vary in the range up to 1 TeV and the_constraint from 
the flavor changing neutral current process & — > 57 is imposeco. The analysis 
of this figure does not include the effects of the spin dependent contribution 
which could change the scatter plot. Specifically, the minimum cross sections 
are sensitive to the inclusion of thc_spin dependent part. Similar analyses can 
be found in other recent reference£3E3. 
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4 Effects of Non-universalities on Dark Matter 

mSUGRA is based on the assumption of a flat Kahler potential. However, the 
nature of physics at the Planck scale is not fully understood. Thus in general 
one should allow for the possibility of a curved Kahler potential. This would 
lead to non-universalities in the scalar sector. However, there are stringent 
constraints on the types of non-universalities allowed in the scalar sector due 
to the FCNC constraint. For example, the FCNC constraint in very strong on 
the amount of non-universality allowed in the first vs the second generation 
sector. However, this constraint is not so strong for the Higgs sector and for 
the third geii£ijation sector. Effects of these constraints have been studied in 
detail in RefsE!l'c3. One finds that in general the presence of non-universalities 
can increase the cross-sections by a factor of 10 or more. 

In addition to modifying the Kahler potential Planck scale physics can also 
modify the gauge kinetic energy function. In general the gauge kinetic energy 
function transforms as the symmetric product of two adjoint representations. 
For the case of SU{5) one finds that the gauge kinetic energy function /^^ 
transforms as the symmetric product of 24 x 24 in SU(5) which contains the 
following irreducible representations of SU{5): 

(24 X 24),y„„„ = 1 + 24 + 75 + 200 (4) 

The SU(5) singlet in the product on the right hand side of Eq.(4) leads to a 
universal gaugino mass while, the additional terms generate non-universalities 
in the gaugino masses at the GUT scale. Thus the SU{3)c x SU{2)l x U{1) 
gaugino masses at the GUT scale are in general admixtures of all the allowed 
representations. This admixture leads to the following relation for the gaugino 
masses at the GUT scale 

mz(0) =mi(l + ^cX) (5) 

r 

where n[ depend on one of the representations on the right hand side in 
Eq.(4) and on the subgroup i. In addition to the appearance of gaugino mass 
non-universalities of the above type in supergravity models, one also finds 
quite naturally non-universalities in a broad class of string models: heterotic, 
Horava-Witten and as well in brane models based on Typel/Type IIB string 
compactifications. Now the value of fi is in general very sensitively dependent 
on non-universalities. For the case of non-universalities in the scalar sector one 
can exhibit in an analytic fashion the dependence of /i on non-universalities in 
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Figure 3: Exhibition of the effects of non-universalities on Ox—v ^ ^ function of the neu- 
traUno mass for different values of C2oo- 



the Higgs sector and in the third generation sector. The analysis shows that /i 
is sensitively dependent on the non-universalities and their effects can signifi- 
cantly decrease the value of and thus affect gaugino vs higgsino composition 
of the neutralino. An_analysis of the effects of non-universalities in the scalar 
sector is given in RefcZI. A similar phenomenon occurs for the case of gaugino 
sector non-universalities. U££e also one can exhibit the dependence of /i on 
non-universalities. One find£3 

where d[i\Jdc2^ > Ojd^j^/dc^^ > 0, 9^200/^^200 < 0. Again with an appro- 
priate choice of the sign of one finds that /i becomes smaller relative to its 
universal value. An analysis o£,the effects of non-universalities in the gaugino 
sector on cr^_p is given in Refts. Fig. 3 gives a typical illustration of the effects 
of gaugino non-universality showing that a significant enhancement of cr^-p 
can occur in the presence of non-universalities. 

5 Bino, Higgsino, and Wino Dark Matter 

The neutralino is in general a mixture of gauginos and higgsinos, 

Xi = XiiB + X12W3 + Xi^Hi + Xi4^H2 (7) 

where B is the bino, W the neutral wino, and Hi,H2 are higgsinos. Now 
an analysis of the neutralino mass matrix shows that for large values of the 
Higgs mixing parameter /i one finds that the neutralino is essentially a bino 
and this situation is realized over a large part of the SUGRA parameter 



6 



spaced. Thus SUGRA models in general predict a bino like cold dark mat- 
ter at least over a major portion of the parameter space. However, in cer- 
tain limited regions of the SUGRA parameter space (i can become small. In 
this case the higgsino and the wino components can become large and one 
may have cold dark matter which is higgsino like or wino like. This phe- 
nomenon can be easily understood by examining the >> Mz limit of 
the components Xi„. Here one finds that the bino component in the neu- 
tralino is given by Xn ~ 1 — (M|/2/i^)sm^0vF, the wino component is of size 
Xi2 ~ (M"^ /2m'^_^p)sin29w sinj3, and the higgsino components have the sizes 
~ —{Mz/fJ,)sin29wsin(3, and X14 ~ {Mz / lj)sin29wsin(3. The above 
illustrates the sensitive dependence of the bino, the wino and the higgsino 
components on /i. One consequence of the effect of large higgsino components 
is that the scalar cross section which depends on the product of the gaugpp 
and the higgsino components increases as the higgsino components increased. 
An analysis ofpCDM with large higgsino components in the context of MSSM 
is given in RefE3 and with large higgsino and wino comnojients in the context 
of anomaly mediated breaking of supersymmetry in RefEU 

6 Comparison of Direct and Indirect Detection 

Indirect detection is complementary to the direct detection in the search for 
dark matter. The most interesting indirect signal for dark matter arises from 
the capture and the subsequent annihilation of the neutralinos in the center 
of the Sun and the Earth. Some of the remnants in the annihilation of the 
neutralinos are the neutrinos which propagate and undergo charged current 
interactions in the rock surrounding the detector and produce upward moving 
muons. The outgoing muon flux can be written in the form — TaI where 
Fa is the xi ~ Xi annihilation rate in the center of the Earth or the Sun and 
f is the product of remaining factors. It, is Ta which is the quantity sensitive 
to SUSY. One can parameterize Ta bjE3 

Ta = ^tanh\t/T) (8) 

Here C is the neutralino capture rate, t is the lifetime of the Earth or the Sun, 
and T = (CC^)~^/^ where Ca is determined by the wimp annihilation cross 
section. An equilibrium between capture and annihilation is reached for the 
case t >> T and in this case one has Ta ^ The equilibrium condition is 
strongly dependent on the susy parameter space. For the case of the sun the 
equilibrium condition is satisfied essentially over all of the parameter space 
while for the Earth it is satisfied only over part of the parameter space. This 
disparity shows up very strongly in the profile of the muon flux from the Earth 
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Figure 4: A comparison of the direct and the indirect detection with a plot of muon flux of 
the Earth and the Sun vs the direct detection rates. Taken from Ref.^ 



and the Sun. In Fig. 4 we exhibit the muon flux for the Earth and the Sun plot- 
ted against the direct detection rate for the germanium detector. The analysis 
of Fig. 4 shows that for relatively large direct detection rates, it is the indirect 
detection from the Earth which is competitive with the direct detection while 
for relatively low direct detection rates it is the indirect detection from the 
Sun which is competitive with the direct detection. Thus indirect detections 
from the Earth and the SuiL^e-complementary. For other recent analyses of 
indirect detection see RefsHaO. 



7 Annual Modulation Effect in Direct Detection 

The annual modulation effect in the direct detection of dark matter is a po- 
tentially important signal for the observation of WIMP like dark matter. The 
effect arises due to the periodicity of the velocity of the Earth ve relative to 
the galaxy,i.e., ve = vs + vocosjcostoit — t^) where vo{= 30fcm/s) is Earth's 
orbital velocity around the Sun, vs{= 232km/ s) is the Sun's velocity rela- 
tive to the galaxy, 7(~ GO'') is the inclination of Earth's orbit relative to the 
galactic plane, and to is June 2. The effects of this motion can produce a 
modulation effect of about 7% in the scattering event rates for wimps. The 
DAMA experiment claims to see an effeciJl. Tlie .- annu al modulation signal has 
been analyzed theoretically in several paperfi3'BEZH2a. In the future one ex- 
pects further data on this exciting possibility. Further, one needs to reduce 
the current ambiguities in the theoretical predictions of dark matter such as 
those discussed in Sec. 2. This requires more accurate chiral perturbation and 
lattice gauge analyses. 
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8 Effects of Large CP Phases on Dark Matter 

It is well known that supersymmetry brings in new sources of CP violation 
since the soft SUSY parameters are in general complex. These CP phases 
induce additional corrections to the electric dipole moments of the electron 
and of the neutron. We consider here the case where the CP phases are large 
and the edm constraiats are satisfied. (For an abbreviated set of references on 
large phases see RefL3). Large phases affect the analysis of dark matter and 
the detailed analyses of event rates for direct detection show that the effeclsjaf 
the phases can change the event rates by an order of magnitude or morcEjO. 
However, with the inclusion of the emd constraints one finds. that the effects 
are significantly reduced although they are still substantiacd. Large phaspi 
also induce a mixing between, the CP even and the CP odd higgs statecBc^ 
which can affect dark mattecJ. 

9 Conclusion 

In this review we have presented the very recent developments in the theo- 
retical analyses of supersymmetric dark matter. These include the effects of 
uncertainties on the wimp density and on the wimp velocity for the Milky Way 
wimps, effects of uncertainties in quark densities in the theoretical predictions 
of cr^-P' ^^'^ effects of non-universalities in the gaugino sector. A anal- 
ysis of the maximum and of the minimum cross sections for x ~ P scattering 
was given. We have also given a comparison of the direct and of the indirect 
detection of dark matter and discussed briefly other recent developments for 
the detection of dark matter. An important topic not discussed is the subject 
of coannihilatioio which can significantly extend the domain of the allasszed 
neutralino masses. In the future one expects that more sensitive detector££lcj 
will be able to probe more deeply the parameter space of SUGRA and other 
SUSY models. 
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